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ABSTRACT 

_ The  flight  enviromnent  of  the  Standard  Missile,  Block  IV  A  interceptor 

involves  high  speeds  that  place  severe  aerodynamic  and  aero-thermal  loads  on  die 
missile  optical  window  and  sensor.  The  supersonic  flow  over  the  missile  produces 
aero-optical  effects  that  can  degrade  IR  seeker  performance.  This  report  summarizes 
an  analysis  of  aero-optical  aberrations  due  to  shock  front  and  shear  layer  density 
gradients.  Computational  Fluid  Dynamics  (CFD)  models,  which  incorporate 
turbulence  due  to  coolant  spray,  are  used  to  compute  the  density  and  refractive  index 
distribution  about  the  IR  dome.  This  data  is  then  used  in  a  ray-tracing  program  to 
determine  boresite  error  and  astigmatic  lensing.  These  results  compare  favorably  with 
computational  results  derived  from  general  fluid  mechanical  approximations. 
Conclusions  are  drawn  which  indicate  that  boresite  error  may  hinder  seeker-targeting 
performance  when  looking  in  a  forward  direction  or  for  low  altitude  intercepts. 


1.0  INTRODUCTION 


The  performance  of  optical  systems  located  in  aerodynamic  environments  can  be 
severelyTmpacted.  These  aero-optical  effects  are  caused  by  index-of-refraction  variations  induced 
by  the  optical  platform  moving  through  the  aerodynamic  flow  field.  Index-of-refraction  variations 
cause  an  optical  system  (both  in  transmission  and  reception)  to  experience  beam  wandCT,  and  beam 
spread  that  results  in  a  r^uction  of  peak  intensity  and  resolution.  These  effects  should  be  considered 
in  the  design  of  airborne  optical  systems  and  performance  evaluation  of  missile  systems. 

Optical  aberrations  resulting  from  aero-dynamic  flow  fields  can  be  categorized  as,  1) 
boresite^ror  and  astigmatic  lensing  due  to  shear  layer  laminar  flow,  and  shear  layer  density  gradient 
about  the  optical  window,  2)  turbulence  in  the  boundaiy  layer,  3)  aero-thermal  effects  resulting  from 
temperature  gradients  across  and  through  the  optical  window. 

For  imaging  in  the  infrared,  a  material  window  suffers  optical  transmission  losses  due  to 
windowradiative  effects  that  induce  background  thermal  noise  on  the  detector  as  well  as  index-of- 
refraction  variations  due  to  temperature  induce  stresses  within  the  window  material.  Aero-thermal 
effects  are  not  considered  in  this  report. 

Flow  testing  data  as  well  as  window  thermal  analysis  support  the  need  for  actively  cooling 
the  IR  window.  Further,  cooling  is  indispensable  for  semiconductor  detector  materials  in  order  to 
maintain  in-band  transmission  and  reduce  thermal  background  emissions.  Conceptually  active 
cooling  maintains  the  window  at  ambient  temperature  preserving  mechanical  properties  such  as 
fracture  strength  and  hardness  of  window  materials. 
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__  There  is,  however,  a  system  cost  to  active  cooling.  Uniform  cooling  of  the  optical  dome 
is  extrenidy  difficult  to  achieve  in  practice.  Computational  Fluid  Dynamics  (CFD)  results  indicate 

_ that  even  through  the  effect  of  the  coolant  gas  is  to  reduce  the  overall  temperature  of  the  dome, 

large  thermal  gradients  develop  across  the  dome  surface,  jeopardizing  material  integrity  and  inducing 
refractive  index  variations  and  hence  optical  aberrations.  Further,  the  coolant  flow,  located  in  the 
dynamic  flow  field,  becomes  turbulent  and  the  refractive  index  variations  in  the  flow  region  result  in- 
pseudo-random  aero-optical  aberrations. 

The  goal  of  this  report  is  to  analyze  optical  aberrations  to  the  missile  seeker,  which  are 
relatively  constant  in  natural  for  a  specific  flight  profile.  These  include  shock  and  shear  layer  effects. 
Aero-optical  effects  resulting  from  the  turbulent  boundary  layer  and  coolant  flow  are  stochastic  and 
necessitate  a  statistical  approach  to  the  analysis,  to  be  summarized  in  a  separate  report. 

_ In  the  next  section,  a  theoretical  basis  for  analyzing  aero-optical  aberrations  due  to  the 

shock  and  shear  layers  is  presented  by  considering  the  refractive  effects  of  an  optical  wave  incident 
on  layCTS  of  slight  index-of-refraction  variation.  CFD  computations,  which  incorporate  turbulence 
due  to  coolant  spray,  are  used  to  compute  the  density  and  refractive  index  distribution  about  the  IR 
dome.  Finally  conclusions  are  drawn  concerning  the  severity  of  aero-optical  aberration,  the 
implications  for  missile  targeting  performance  and  suggestions  for  further  work. 


2.0.  THEORY 


2.1.  SHOCKWAVES 

_ Shock  waves  result  as  a  body  accelerates  through  the  atmosphere  from  the  subsonic  into 

the  supersonic  region.  At  supersonic  velocities  a  shock  wave  develops,  as  discontinuities  in  the  air 
properties  such  as  density. 

Due  to  the  density  discontinuity  and  hence  index-of-refraction  variations,  the  shock  front 
both  refracts  and  reflects  the  incoming  optical  radiation  depending  on  the  strength  and  angle  of  the 
shock  front.  In  addition  shock  fronts  usually  cause  low-order  aberrations.  For  the  case  of  a  rapidly 
maneuvering  body,  the  aero-optical  induced  aberrations  can  vary  with  time,  typically  with  a  temporal 
frequency  related  to  the  duration  of  the  maneuver. 

The  relation  between  the  density  ahead  of  the  shock  front,  pi,  and  density  behind  the 
front,  given  by;' 


P2  _  (7 -1- 1)A/o  sin^  0 

Pi  ~  (7-l)Mosin^0  +  2 


where  y  is  adiabatic  index 

Mo  is  free-stream  Mach  number 
0  is  the  shock  angle  at  the  optical  axis 

j  No  analytical  scaling  relationships  have  been  developed  for  determining  shock  front 
angle^Fdensity  gradient;  however,  approximations  have  been  derived  based  on  CFD  computations 
The  strength  of  the  shock  as  measured  by  the  density  gradient  across  the  shock  increases  with 


*  Landau  and  Lifshits,  Fluid  Mechanics,  Vol.  6,  pg.  331,  Pergamon  Press,  1959. 

^  National  Advisory  Committee  for  Aeronautics  (NACA),  Washington,  DC,  “Equations,  Tables,  and 
Charts  for  Compressible  Flow”  (U).  Technical  Report  1135,  1947. 


Mach  number.  The  angle  of  the  shock  front  depends  on  the  angle  of  attack,  geometry  of  the 
vehicle,  as  well  as  the  Mach  number.  Blunt  bodies  have  shock  detached  and  standing  ahead  of  the 
body.  Surfaces  with  sharp  leading  edges  exhibit  shock  waves  attached  to  the  leading  edge.  Figure 
shows  representative  curves  of  shock  wave  angle  as  a  function  of  the  cone  angle  of  die  body. 


Figure  1.  Shock  Wave  Angles. 


2.2.  TRACKING  OR  BORESITE  ERROR 

Density  values  p,  which  result  from  aero-dynamic  analysis  are  related  to  refractive  index 
values,  nThy  the  Gladstone-Dale  law,^  n  - 1  =  Gp,  where  G  is  the  Gladstone-Dale  parameter.  The 
Gladstone-Dale  parameter  is  plotted  in  Figure  2  and  can  be  seen  to  be  nearly  constant  at  0.22  cmVg 
for  the  wavelength  X  >  0.5  pm.'' 

A  complex  shock  and  flow  pattern  that  is  an  unique  function  of  the  missile  speed, 
attitude,~^tude,  and  coolant  mass  flow  rate  is  formed  upstream  of  the  dome.  A  shock  wave  is 
usually  very  stable  for  constant  flight  characteristics.  Since  the  shock  front  is  usually  at  some  angle 
to  the  optical  axis,  an  incoming  light  ray  is  refracted  at  the  shock  front  interface.  Using  SnelTs  Law 
and  the  Gladstone-Dale  relation,  the  angle  of  the  incident  ray.  Pi,  is  related  to,  P2.  the  angle  of  the 
refracted  ray,  both  measured  normal  to  the  shock  plane,  by^ 


*  Philosophical  transactions  of  the  Royal  Society  of  London.  Series  A,  Containing  papers  of  a 
mathematical  or  physical  character.  London,  Royal  Society  of  London,  Vol.  148,  pp.  882  or  887. 

G.  W.  Gladstone,  et  al.,  “On  the  Influence  of  Temperature  on  the  Reflectivity  of  Light  (U).” 

“*  Atmospheric  Propagation  of  Light,  The  Infrared  and  Electro-Optical  Systems  Handbook  (U),  Vol.  2, 
pg.  254,  ERIN,  SPIE  Press. 

’William  L.  Wolfe  and  George  J.  Zissis,  Eds.  The  Infrared  Handbook  (U),  The  Infrared  Information 
and  Analysis  (IRIA)  Center,  Environmental  Research  Institute  of  Michigan,  1978. 


(2) 


sin  ^  =  sin  Pi 


l  +  Cpl' 

1  +  <^P2  y 


Figure  2.  _ The  Gladstone-Dale  Parameter  as  a  Function  of  Wavelength. 

Note  that  for  X  >  0.5  |jm,  this  parameter  is  nearly  invariant  (0.22  cmVg). 


2.3.  FOCAL  LENGTH  VARIATIONS 

_ The  shock  front  can  be  represented  by  laminar  density  contours  in  the  plane  of  incidence 

formed  by  the  optical  ray  and  missile  axis.  In  the  orthogonal  direction  the  shock  front  is  curved  and 
acts  as  an  aerodynamic  cylindrical  lens,  producing  astigmatism  in  that  the  sagittal  ray  focus  behind 
the  tangential  rays.  Figure  3  illustrates  this  case.  The  curved  shock  acts  like  a  lens  of  radius  equal  to 
the  distance  from  the  vehicle  center  to  the  shock  front.  The  effective  focal  length  F  in  the  plane 
perpendicular  to  the  airframe  longitudinal  axis  and  at  station  distance  x  is,® 


F  _  I  Ml  sin^  0(1  +  6Gpi )  +  5 
X  6Gp]Mi  sin^  6 


As  can  be  seen  from  the  above  equation,  the  lensing  effect  becomes  stronger  with 
increasmg^ach  number  and  cone  angle.  The  effect  of  the  conical  lens  is  to  skew  the  focal  plane, 
ran  sing  a  defocus  of  the  optical  spot.  The  lens  effect  is  actually  a  more  complicated  mechanism  in 
that  the  radius  of  the  shock  depends  upon  the  sensor  look  angle  since  the  income  optical  ray  takes  a 
different  path  from  the  shock  front  to  the  sensor.  When  considered  in  this  manor,  the  shock  front  is 
parabolic  in  sh^e  and  requires  a  complex  integration  for  complete  analysis.  Fortunately,  the 
focusing  shift  due  to  the  curved  nature  of  the  shock  front  is  small  and  the  spherical  shock  surface 
approximation  is  adequate. 


*  Atmospheric  Propagation  of  Light,  The  Infrared  and  Electro-Optical  Systems  Handbook  (U),  Vol, 
2,  pg.  252,  ERIN,  SPIE  Press. 


Figure  3  Shock  Wave  Induced  Lensing  Effect. 


If  we  consider  the  aero-dynamic  lensing  of  the  shock  front  and  the  optical  imaging 
system  ^wo  thin  lenses  then  the  effective  focal  length  of  the  system  F.  is  related  to  the  shock  front 
focal  length  F  and  system  focal  length  F„,  by 


1  1 
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2.3.  SHEAR  LAYER  OPTICAL  EFFECTS 

The  shear  layer  which  forms  from  the  aerodynamic  flow  about  the  optical  dome  also 
contribute  to  refraction  and  lensing  of  incident  optical  radiation.  CFD  computations  reveal  a  density 
gradient  across  the  shear  layer  typically  greata:  than  that  of  the  shock  front.  Further,  the  shear  layer 
is  curved  in  the  plane  of  incidence  as  the  flow  conforms  to  the  hemispherical  dome,  producing  a 
lensing  and  focal  shift  in  the  plane  of  incidence  as  well  as  the  orthogonal  plane. 

Tlie  influence  of  the  shear  layer  to  the  optical  system  is  considered  by  two  separate 
technique  First,  a  computational  technique,  similar  to  that  used  in  analyzing  boresite  error  due  to 
the  shock  front,  is  extended  to  include  the  shear  layer.  In  this  method,  the  incoming  optical  beam  is 
first  refracted  by  the  shock  front  and  intersects  the  shear  layer.  The  slope  and  strength  (density 
gradient)  of  the  shear  layer  is  determined  by  CFD  results  and  depends  on  the  optical  sensor  look 
angle.  Snell’s  law  (eg.  2)  is  then  used  to  find  the  refractive  angle  through  the  shear  layer  and  the 
overall  boresite  error  due  to  refractive  effects  from  the  shock  and  shear  layer  determined. 

In  the  second  method  optical  effects  induced  by  the  shear  layar  and  shock  front  are 
analyzed  by  employing  CFD  density  data.  Hare  the  density  contours  are  modeled  as  thin  wedges  of 
varying  optical  thickness  and  a  ray  tracing  code  used  to  determine  the  magnitude  of  various  optical 
distortions.  This  technique  is  described  in  more  detail  in  the  following  section. 


2.4.  RAY-TRACING  FROM  COMPUTATIONAL  FLUID  DYNAMICS 


In  this  section,  CFD  models,  which  inco^orate  turbulence  due  to  coolant  spray,  are  used 
to  compule  the  density  and  refractive  index  distribution  about  the  IR  dome  including  the  shock, 
boundary,  and  shear  layers.  CFD  density  plot’  for  a  nominal  flight  profile  is  shown  as  Figure  4.  This 
data  is  then  used  in  a  ray-tracing  program  to  determine  boresite  error  and  astigmatic  lensing.  These 
results  compare  favorably  with  computational  results  derived  from  general  fluid  mechanical 
approximations  reported  above. 


Figure  4.  CFD  Generated  Density  Contours  in  0.05  kg/m^  Steps  for  Nominal  Fight  Profile, 


Ray-tracing  from  CFD  results  are  briefly  described  for  a  specific  look  direction.  The 
CFD  denHty  contours  in  steps  of  0.05  kg/m^  are  converted  into  index-of-refraction  (n)  contours  via 
the  Gladstone-Dale  relation.  From  the  slope  of  the  density  contour  lines  with  respect  to  the  look 
angle,  and  the  thickness,  d,  of  each  layer  between  contours,  the  air  is  modeled  as  thin  wedges  of 
appropriate  optical  thickness,  nd,  and  tilt  angle.  The  layers  were  given  curvature  in  the  orthogonal 
plane  to  model  the  lensing  effect.  The  CFD  density  contours  modeled  as  an  optical  system  are 
illustrated  in  Figure  5.  Index  of  refraction  and  wedge  angle  data  are  used  to  model  the  flow  field  in 
an  optical  ray-tracing  code  (ZEMAX).  A  2.286-cm  diameter  entrance  stop  was  placed  on  the 
window  rather  than  at  the  sensor  entrance  pupil.  A  perfect  lens  of  5. 75 -cm  focal  length,  matching  the 
optical  system  of  the  seeker,  was  used  in  the  ray-tracing  routine.  Three  look  angles  were  used,  0° 
(perpendicular  to  the  missile  longitudinal  axis),  45°  and  80°. 

Results  of  the  geometric  ray-trace  are  shown  as  boresite  error  and  focal  spot  shape, 
althouglf^any  other  performance  evaluation  tools  are  available.  These  techniques  most  clearly 
indicate  optical  degradation  due  to  the  shock  and  shear  layer. 


Ron  Schultz,  Naval  Air  Warfare  Center  Weapons  Division,  China  Lake,  Calif. 


ZERO  DEGREE  LOOK  ANGLE/NOMINAL  FLIGHT  PROFILE 


Plane  of  Incidence 


Orthogonal  to  Plane  of  Incidence 


Figure  5  CFD  Density  Profiles  Modeled  as  Optical  Wedges  of  Varying  Index. 

3.  RESULTS 

The  results  of  this  paper  are  divided  into  the  two  sections.  First,  results  derived  from  the 
general  shock  front  theory  are  presented.  This  includes 

1  Boresite  error  as  a  function  of  look  angle  for  nominal  flight  condition. 

2.  Boresite  error  as  a  function  of  missile  altitude. 

3.  Aerodynamic  flow  induced  lensing, 

Ray-tracing  analysis  from  CFD  generated  data  are  shown  next.  This  analysis  includes, 

1  Boresite  error  as  a  function  of  look  angle  for  nominal  flight  conditions. 

2.  Lensing  induced  astigmatism. 


To  detamine  boresite  error.  Equation  (1)  was  first  used  to  calculate  the  density  gradient  across 
die  shock  front.  A  missile  corie  angle  of  20°  was  used  in  Figure  1  to  arrive  at  a  17°  shock  wave  angle, 
which  compares  favorably  to  CFD  results.  Diffractive  effects  due  to  the  shear  layer  were  based  upon 
density  gradients  determined  from  CFD. 

Boresite  error  results  are  shown  in  Figures  6  and  7. 
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Figure  6. _ Boresite  as  a  Function  of  Missile  Attitude 

for  Nominal  Flight  Profile. 


Figure  7.  _ Boresite  for  Nominal  Flight 

Profile  as  a  Function  of  Look  Angle  Zero 
Degree.  Angle  is  straight  up.  90  degrees  is 
looking  forward. 


_ Figure  6  plots  boresite  error  for  a  missile  traveling  at  Mach  4.1  at  10-degree  angle  of 

attack.  The  sensor  look  angle  is  45  degrees.  Boresite  error  is  inversely  proportional  to  the 
atmospheric  density,  since  the  density  gradient  across  a  shock  front  is  directly  proportional  to  the  free 
stream  density  gradient. 

_ Boresite  error  as  a  function  of  look  angle  for  a  nominal  flight  profile  of  15-km  altitude, 

Mach  4.1  and  10-degree  angle  of  attack  is  shown  in  Figure  7.  Results  for  both  the  general  theory 
and  ray  tracing  from  CFD  results  are  shown.  The  boresite  error  is  quite  small  at  small  look  angles  but 
rises  rapidly  above  45  degrees.  The  focus  shift  for  a  seeker  under  nominal  flight  conditions  was 
determined  from  Equations  3  and  4  and  plotted  in  Figure  8  as  a  function  of  look  angle.  This 
calculation  assumes  the  optical  system  is  focusing  on  a  target  10-km  distance.  As  the  look  angle 


approaches  40  degrees  the  focus  shift  effect  increases  rapidly,  but  the  overall  effect  is  still  small 
^iice  the  astigmatic  spot  size  is  well  undo:  the  diffraction  limited  spot. 


Figure  8.  Effective  Focal  Distance  as  a 
Function  dfXook  Angle  for  a  Nominal  Fhght 
Profile. 


Boresite  oror  determined  from  ray  tracing  based  on  CFD  data  is  also  plotted  in  Figure  7 
at  three  Took  angles.  Good  agreement  between  the  general  theoretical  predictions  and  ray  tracing  was 
obtained  for  small  look  angles.  Ray  tracing  resulted  in  smaller  boresite  error  than  predicted  by  the 
general  theory  as  the  look  angle  increased. 

Spot  diagrams,  generated  by  ray  tracing,  are  shown  in  Figures  9,  10,  and  11  for  0, 45,  and 
85-degre^ook  angles,  respectively.  These  diagrams  result  from  tracing  numerous  rays  through  the 
aero-optical  system,  neglecting  diffraction  effects,  and  displaying  the  pattern  that  would  result  on  the 
sensor  detector  plane.  As  can  be  seen,  a  large  amount  of  astigmatism  is  present  in  the  spot  diagrams, 
but  the  extent  of  the  spot  elongation  is  relatively  small;  0.4,  1.0,  and  2.0  pm  for  zero,  45  and  85 
degree  look  angles,  respectively.  These  values  of  lensing  induced  astigmatism  are  small  compared  to 

the  Airy  diffi-action  spot  of _ pm  representing  the  diffraction  limit  of  the  optical  system  with  0.9- 

inch  entrance  diameter  and  2.25  inch  focal  length.  The  final  diffraction  pattern  at  the  optical  systems 
focal  plane  will  be  a  convolution  of  the  astigmatic  aberrated  spot  and  the  diffraction  limited  airy 
pattern. 
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Figure  9. 


Zero  Degree  Look  Angle 
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Figure  10. 


45  Degree  Look  Angle. 
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Figure  II  85  Degrees  Look  Angle. 

3.1  OPTICAL  REFLECTION  LOSS 

Light  interacting  with  the  laminar  shock  front  is  both  refracted  on  transmission  and 
reflected.  The  reflected  component,  R,  for  unpolarized  Ught  is  given  by 


From  the  Gladstone-Dale  relations  and  assuming  no  =  1,  the  reflectance,  R,  is; 


'  Gp  f 

^Gp  +  lj 


The  reflectance  is  on  the  order  for  10  ®  for  all  look  angles  and  is  therefore  considaed 

negligible 


4.  CONCLUSIONS 


_ This  report  summarizes  the  analysis  of  aerodynamic  induced  optical  distortions  to  the 

Standard  Missile  Block  IV  A  interceptor. 

Interaction  of  the  flow  field  with  the  airborne  platform  and  sensor  window  can  induce  a 
variety  ofoptical  effects.  Density  gradients  in  the  shock  front  and  shear  layer  about  the 
hemispherical  window,  boundary  layer  turbulence,  window  emissivity  and  temperature/stress  gradients 
result  in  optical  abarations,  which  potentially  jeopardize  missile  tracking  and  targeting  performance. 
This  report  considers  aberrations  resulting  from  the  shock  and  shear  layer  density  only.  Optical 
effects  do  to  random  aerodynamic  turbulence  and  aero-thermal  induced  aberrations  are  reported 
separately. 


The  analysis  undertaken  in  support  of  this  project  was  two-fold,  first  a  theoretical 
approaclTbased  upon  general  fluid  mechanical  theory  and  second,  CFD  models,  which  incorporate 
turbulence  due  to  coolant  spray  used  in  a  ray  tracing  program  to  determine  boresite  error  and 
astigmatic  lensing.  These  results  compare  favorably  with  computational  results  derived  from  general 
fluid  mechanical  approximations. 

Specific  conclusions  follow: 

Line-of-sight  (boresite  orors)  are  small  for  a  Mach  4.1  interceptor  at  15  km  altitude 
afid  zero  degree  (up)  sensor  look  angle.  This  is  equivalent  to  30-cm  error  at  10-km  target 
range. 

2.  As  the  look-angle  increases,  however,  boresite  becomes  more  severe,  approaching  2 
meters  at  10-km  distance  and  78°  look-angle. 

3.  Boresite  error  also  decreases  with  increasing  altitude,  as  expected.  This  has 
i^lications  for  interception  at  low  altitudes.  Boresite  error  for  a  target  10-km  distance 
and  2.5-km  intercept  altitude  is  estimated  at  4  meters. 

4.  _ Lensing  induced  astigmatism  is  negligible,  generally  much  smaller  than  the  diffraction 

limited  spot. 

5.  Reflection  loss  at  the  shock  interface  is  negligible. 


